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We emphasise the importance of the gas entropy in studying the evolution of cluster 
gas evolving under the influence of radiative cooling. On this basis, we develop an 
analytical model for this evolution. We then show that the assumptions needed for 
such a model are consistent with a numerical solution of the same equations. We 
postulate that the passive cooling phase ends when the central gas temperature falls 
to very low values. It follows a phase during which an unspecified mechanism heats the 
cluster gas. We show that in such a scenario the small number of clusters containing 
gas with temperatures below about 1 keV is simply a consequence of the radiative 
cooling. 



1. Introduction 

Various speakers at this meeting have shown over- 
whelming evidence that to explain our observations of 
the hot gas in galaxy clusters we need to invoke some 
form of non-gravitational heating of this gas. Various 
mechanisms for heating the gas in galaxy clusters have 
been proposed. All of these can be divide into two cat- 
egories: Steady-state or episodic heating scenarios. In 
the steady-state picture the heating mechanism 'knows' 
locally about the radiative cooling. Every bit of energy 
radiated away by the gas at a given place inside the clus- 
ter is replaced by an equal bit of energy at the same 
place by the heating process. The cluster as whole never 
changes. For episodic heating this tight spatial and tem- 
poral connection between cooling and heating is broken. 
The cluster gas spends some time passively cooling and 
the lost energy is then replenished during a heating pe- 
riod. The energy input does not have to be distributed in 
this case. For example, the heating of the cluster gas by 
an AGN outflow may initially only affect the very centre 
of the cluster. However, such a localised release of energy 
adds the further complication of having to subsequently 
distribute the energy throughout at least the inner vol- 
ume of the cluster. This distribution of energy cannot be 
achieved by a fully convective flow as this would imply a 
negative entropy gradient in the cluster which has never 
been observed. 

In this contribution we are not investigating any heat- 
ing mechanisms. We simply assume that the cluster gas 
is heated in an episodic fashion, possibly by AGN out- 
flows. In between heating phases the cluster gas cools 
passively until the temperature of the gas at the cluster 
centre vanishes. At this point the heating mechanism is 
triggered once more. In the following we solve the hy- 
drodynamical equations governing the evolution of the 
cluster gas during the passive cooling periods. We will 
show that the properties of the cluster gas predicted by 
our model are consistent with observations. In particular 
we demonstrate that our model can explain why we find 
very few clusters containing gas at temperatures below 



about IkeV. 

In section|21we motivate our approach of concentrating 
on the entropy of the cluster gas rather than its internal 
energy. Section |31 outlines the development of a simple 
analytical model describing the evolution of the cluster 
gas. The assumptions for the analytical model are tested 
against a numerical solution of the same equations in 
section 0] Finally, in section [S] we use the numerical 
approach to demonstrate that passive radiative cooling 
alone can explain the apparent paucity of detections of 
clusters containing very cold gas. 

2. Why is entropy important? 

In the following we give two reasons why the entropy 
of a gas is an important quantity to study in galaxy clus- 
ters. The first is a significant simplification of the hydro- 
dynamic equations describing the cooling flow gas. Sec- 
ondly, observations of cluster gas reveal a simple power- 
law relation between the gas' entropy and its mass. 

2.1. Comparison with energy 

The properties of the hot gas in galaxy clusters evolv- 
ing under the influence of radiative cooling are tradition- 
ally described by a set of three hydrodynamical equa- 
tions: Conser vation of mass, momentum and energy (e.g. 
ISarazinlfl98i^ . Usually these equations are then solved 
numerically and the resulting inflow of gas towards the 
cluster centre due to radiative cooling is referred to as a 
cooling flow. It is intuitive to use the conservation of en- 
ergy equation as the rate of change of the internal energy 
of the gas due to radiative coohng is very simple, i.e. 



dE, 



rad 



dt 



(1) 



where p is the mass density and A(r) is the temperature- 
dependent cooling function. However, the total rate of 
change of the internal energy is much more complicated. 
As the gas cools, it is compressed by the somewhat hotter 
gas further out in the cluster. This adiabatic work in- 
creases the internal energy. Radiative cooling is a strong 
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function of the gas density. Thus coohng is most impor- 
tant in the cluster centre and the continued compression 
of the gas there leads to a gas flow inwards. On the way 
in the gas looses gravitational potential energy. Other 
ways in which the gas may gain or loose energy include 
the conversion of kinetic energy of the gas flow and ther- 
mal conduction. All these processes make the equation 
of energy difficult to solve. 

In fact, we can formulate an equivalent, but consider- 
ably simpler set of equations to describe the evolution 
of the cluster gas. For this, we consider a variant of the 
entropy of the gas defined as 

S = Nks const. + ^ In (pp 

where N is the number of gas atoms, fee is the Boltz- 
mann constant and p is the gas pressure. Here we have 
assumed a monatomic gas with adiabatic index equal to 
5/3. Entropy itself is somewhat cumbersome to work 
with and so we define the 'entropy index' as 

Entropy and therefore the entropy index does not 
change for adiabatic compression of the gas. If we further 
assume that all gas motions are subsonic, i.e. the kinetic 
energy of the gas is small compared to other forms of 
energy, and that thermal conduction is suppressed, then 
the only process changing the gas entropy is radiation. 
It is straightforward to show that 

d = -^pi/3A(r). (4) 

We will see in the following that using this simple expres- 
sion instead of the conservation of energy equation leads 
to an analytical solution of the hydrostatic equations. 

2.2. Entropy index - mass relationship 

iKaiser fc BinnevI l)2003() pointed out that in the Hydra 
cluster the cumulative mass of gas with an entropy index 
less than a given cr is a simple power-law function of a 
itself, i.e. 

M{<<j)= A {a- <Jo)\ (5) 

where ^ is a constant and ctq is the entropy index of the 
gas with the lowest entropy in the entire atmosphere lo- 
cated at the centre of the cluster. Figure ^ shows how 
closely a power-law of this form fits the observatioi is. An- 
other example is the Virgo cluster (see Figure El iKaised 
More clusters for which a similar relation may 
very well provide an excellent fit to the data can be found 
in the contributions to these proceedings by Horner et al. 
and Donahue. 

A parcel of gas with a given entropy index at a specific 
time t has a well-defined mass. This mass will not change 
under the influence of radiative cooling. Therefore at 
time t + dt the parcel of gas will have a lower entropy 
index, but the mass of the gas in the parcel with this new 
value of cr will not have changed. Thus we can suspect 
that the power-law relation between M{< a) and a — (Tq 
may persist under the influ ence of radiativ e coo ling. This 
was in fact confirmed bv lKaiser fc BinnevI l)2003 ,. see also 
section 0J by numerical integration of the hydrodynami- 
cal equations. An analytical argument for this behaviour 
will be presented in a forthcoming paper. 
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Fig. 1. — Distribution of gas mass as function of entropy index, 
/ , for the Hydra cl uster. Triangles show observational results from 
The line shows the best-fitting relation of the 
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form given by equation |^ . The units of cr are cm^ g 



2/3, 



101" 




cr / keV cn 



Fig. 2. — Distribution of gas mass as function of entropy index, cr, 
for the Virgo clus ter. Data points show observational results from 
IMatsushita et ahl ( 2002). The line shows the best-fitting relation 
of the form given by equation JSJ. 



3. Analytical model for evolving cooling flows 

We use equation Q giving the rate of change of the 
entropy index to replace the conservation of energy equa- 
tion. The conservation of momentum equation remains 
unchanged, i.e. 



dp 
dr 



9$ 
or 



p\3/5 a$ 



(6) 



where $ is the gravitational potential. As usual wc as- 
sume that $ is generated by the dark matter halo of the 
cluster and remains unchanged by the motion of the gas. 

We now assume a cooling function A(T) appropriate 
for radiative cooling due to pure bremsstrahlung. This 
implies A(r) cx yT. Clearly this is a strong simplifica- 
tion as it neglects any cooling due to line emission. Even 
for a gas of vanishing metallicity line emission will change 
the cooling function significantly for gas with tempera- 
tures below about 2 keV. However, above 2 keV this ap- 
proximation of the cooling function is reasonable. Due 
to the very efficient cooling at low temperatures most of 
the cluster gas will spend only very little time with a 
temperature below this threshold. With this we find 



aocp2/V/i". 



(7) 
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Fig. 3. — Density distribution of the cluster gas in the analytical 
model for two different gravitational potentials. 
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Fig. 4. — Temperature distribution of the cluster gas in the 
analytical model for two different gravitational potentials. 



Clearly the rate of entropy loss is only a weak func- 
tion of the changing gas properties. Therefore we set 
p2/5^i/io const. Apart from the factor a-^^^^, this is 
equivalent to the assumption of isobaric cooling roughly 
consistent with the fact that the cooling timescale is long 
compared to the dynamical timescale of the gas. 

Our final assumption in order to facilitate an analytical 
solution is that of a uniform entropy index, cTj, of all the 
gas throughout the cluster atmosphere aX t = 0. Again 
this may well be a gross oversimplification of the real sit- 
uation as the cluster gas is presumably accumulated from 
a number of dark matter subhaloes merging to form the 
final cluster halo. There is no reason why the gas con- 
tained in these subhaloes should have the same entropy. 
However, any heating process affecting the entire clus- 
ter will eventually lead to a uniform entropy throughout 
the cluster atmosphere. Whether such a comprehensive 
heating phase has taken place or not is of course not 
clear. However, the assumption is crucial to obtaining 
an analytical solution. 

With these assumptions we can solve equations Q and 
® to find 

a^{ar^comf' (8) 



and 



P = 



(a, - a) - 



5/2 



(9) 



Here, co and c\ are constants and the other thermody- 
namic properties of th gas may be derived from the above 
two expressions by assuming ideal gas conditions. Also 
note that the solution does not depend on the form of 
the gravitational potential $. 

Figures 13 and ^illustrate the density and temperature 
distributions in the cluster atmosphere predicted by the 
analytical model. As expected, assuming a /3-profile for 
the gravitational potential results in the formation of a 
uniform density core while a NFW-profile ijNavarro et al.l 
[l99 6) results in a steeper central density distribution. 
The temperature distributions show an off-centre peak 
for both distributions, but the /?-profile also results in 
a constant temperature core. This general shape with 
a constant temperature followed by a peak further out 
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Fig. 5. — Cumulative mass, M(< tr), as a function of entropy 
index as predicted by the analytical model. The vertical lines show 
the location of the gas within the cluster for the model with a 
gravitational potential described by the /3-profile. 



is reminiscent of the gener i c ten iperature profile found 
empirically bv lAllen et alJ l)200lD studying a sample of 
cooling flow clusters. 

Finally, Figure |31 illustrates that, over a large range in 
entropy index and radius, the analytical model predicts a 
power-law relation between the gas mass and the entropy 
index as described by equation (O. The slope of the 
power-law, e, in the inner part of the relation is roughly 
equal to 3 for the NFW-profile and 3/2 for the /3-profile 
independent of the value of /3. Given that the relation 
steepens for increasing radius, the latter value is closer 
to the slopes (Hydra: e ^ 1.5, Virgo: e ~ 2.3) found 
empirically from the data shown in Figures ^ and |21 

4. Numerical approach 

The results from the analytical model presented above 
may arise as a direct consequence of the great simplifica- 
tions we employed in deriving this solution. In order to 
test this possibility, we have also solve equations Q and 
10 by numerical integration. Without the assumptions 
we made above, we need a third equation to close the 
system. Here we use a variant of the equation of mass 
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Fig. 6. — Entropy index - mass relation from the numerical 
solution. The long-dashed line shows the relation at the beginning 
of the integration approximating the observational data obtained 
in the Hydra cluster (see Figure 0. The relation at the end of 
the calculation when the central entropy index approaches zero is 
shown by the solid line. The short-dashed line shows the power- 
law relation of the form of equation j^J that provides the best 
approximation to the solid line. 



conservation in the form 



dM _ DM da 
dr da dr 



(10) 



We start the integration with the gas distribution de- 
rived fr om the X-ray obser vations of the Hydra cluster 
given in 'David et alJ l)200lh . These data provide us with 
the initial values for the derivative dM/da. We then use 
equation I^J to evolve the entropy index through a small 
timestep dt. We calculate the new derivative using the 
fact that the mass of gas with a given entropy index does 
not change (see section lT^ . Now it is possible to numer- 
ically solve equations l@J and (|10() . We then repeat the 
whole procedure for the next timestep. The numerical 
solution also includes a more realistic cooling function 
incorporating linc emission. Details of the calculations 
are given in lKaiser fc BinnevI lj200^ . 

We continue the calculation until the entropy index 
and the temperature of the gas at the cluster centre are 
equal to zero. Figure demonstrates that the entropy 
index - mass relationship continues to be well approxi- 
mated by a power-law throughout the calculation. The 
inclusion of radiative cooling due to line emission does 
not change this relationship. 

The entropy evolution is strongest at the cluster centre 
as the pressure is highest there. Therefore we expect the 
deviations of the entropy evolution from a strictly linear 
behaviour to be strongest in the cluster centre. Figure 
[7| shows the evolution of the entropy index of the gas at 
the cluster centre. Except for very late times the graph 
shows that this evolution is linear with time. Thus we 
conclude that our assumption for the analytical model 
of a linear evolution of the entropy index throughout the 
cluster is a good approximation to the real situation. 

5. Why we do not detect cold gas in clusters 

A further prediction of the model in both its analytical 
or numerical form is the strongly accelerating evolution 
of the gas temperature in the cluster centre. Figure |H1 
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Fig. 7. — The evolution of the entropy index of the gas at the 
cluster centre in the numerical solution. 




100 1 50 200 
Time / Myrs 

Fig. 8. — The evolution of the temperature of the 
cluster centre in the numerical solution. 
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shows the evolution of the gas temperature at the clus- 
ter centre. By construction this is the coldest gas in the 
entire cluster. The temperature decrease only gradually 
for a few 10^ years to about 2keV. Only then the tem- 
perature drops rapidly to below 1 keV. This rapid drop 
is enhanced by the increasing importance of cooling due 
to line emission at about this temperature, but the drop 
would also occur in the absence of line cooling. The im- 
portant point is that the lowest gas temperature in the 
cluster exists only for a very short fraction of the clus- 
ter evolution timescale. Furthermore, the volume of gas 
with temperatures below 1 keV or even 2 keV is small. 
In other words, the total emission measure of gas with 
temperatures below a given low threshold is only a small 
fraction of the total emission measure of the cluster. 

Consider a sample of clusters with a given lower flux 
limit. Assume that every cluster in the sample has the 
same properties as the Hydra cluster at i = 0, but is ob- 
served at a random time in its evolution. We can then use 
the numerical model to estimate the fraction / of clus- 
ters in this sample which contain a detectable amount of 
gas at temperatures below a given threshold. Figure |51 
presents the results of such an estimation. Even if the 
original survey is followed up by much more sensitive 
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Fig. 9. — The fraction / of clusters in a flux-liinitcd sample 
that contain a detectable amount of gas with temperatures below 
a threshold T as a function of T. The lines show the fraction / 
for various flux limits used to detect the cold gas in follow up ob- 
servations to the original survey. Dotted line: The same flux limit 
as the original survey. Solid lines from bottom to top; 1000, 2000 
and 10** times lower flux limits compared to the original survey. 



observations to search for cold gas, the fraction of clus- 
ters with detectable amounts of gas with temperatures 
below 1 keV is very small. It is therefore not surprising 
that only very few clusters containing cold gas have been 
found. 

6. Summary 

In this contribution we studied the evolution of the hot 
gas in galaxy clusters under the influence of radiative 
cooling. By concentrating on the entropy of the gas, we 



were able to find an analytical solution for the relevant 
hydrodynamical equations. This solution assumes that 

• the radiative cooling is due to pure bremsstrahlung. 

• the evolution of the entropy index of a given parcel 
of gas is linear in time. 

• the entropy index has a uniform value throughout 

the cluster atmosphere at t = 0. 

With the help of a numerical integration of the same 
basic equations we showed that all of these assumptions 
are reasonable. 

Both the analytical model and the numerical solution 
predict a simple power-law relationship between the en- 
tropy index and the gas mass consistent with observa- 
tions. The accelerating temperature evolution of the 
model can explain why we do not observe many clusters 
with very cold gas at their centres. However, this result 
depends crucially on the triggering of a heating mech- 
anism for the cluster gas at the time when the central 
temperature in the cluster approaches very low values. 
Without such a heating mechanism, this model encoun- 
ters the usual problem that cold gas would rapidly accu- 
mulate at the cluster centre. How the required heating 
mechanism works is not clear. However, the arrival of 
cold gas at the cluster centre at the time when the mech- 
anism must be triggered may point to a close connection 
between these two events. 



CRK would like to thank PPARC for financial support. 



References 



Allen, S. W., Schmidt, R. W., & Fabian, A. C. 2001, MNRAS, 328, 
L37 

David, L. P., Nulsen, P. E. ,1., McNamara, B. R., Forman, W., 
Jones, C, Ponman, T., Robertson, B., & Wise, M. 2001, ApJ, 
557, 546 

Kaiser, C. R. 2003, MNRAS, 343, 1319 

Kaiser, C. R. & Binney, J. J. 2003, MNRAS, 338, 837 



Matsushita, K., Belsole, E., Finoguenov, A., & Bohringer, H. 2002, 

A&A, 386, 77 

Navarro, J. F., Frenk, C. S., & White, S. D. 1996, ApJ, 462, 563 
Saraain, C. L. 1988, X-ray emission from clusters of galaxies 
(Cambridge University Press, Cambridge.) 



